On the universality of the scahng of fluctuations in traflic on complex networks 
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We study the scahng of fluctuations with the mean of traffic in complex networks using a model 
where the arrival and departure of "packets" follow exponential distributions, and the processing 
capability of nodes is either unlimited or finite. The model presents a wide variety of exponents 
between 1/2 and 1 for this scaling, revealing their dependence on the few parameters considered, 
and questioning the existence of universality classes. We also report the experimental scaling of the 
fluctuations in the Internet for the Abilene backbone network. We found scaling exponents between 
0.71 and 0.86 that do not fit with the exponent 1/2 reported in the literature. 

PACS numbers: 



Reeently, the theory of complex networks has started 
to cope with the problem of dynamics on networks. After 
much work devoted to the understanding of the network 
topology 0, the physics community has begun to de- 
velop models to uncover the phenomena observed in the 
dynamics on complex networks. Some stylized models of 
traffic flow in complex networks HHIilEllIllilillS 
mi can be used to gain intuition about complex networks 
dynamics, and to determine the leading parameters of the 
dynamic processes related to the network topology. 

The main results obtained up to now concerning traffic 
flow in complex networks, are related to the determina- 
tion of bounds for this flow to become congested. Never- 
theless, traffic on real complex networks, as for example 
the Internet, is not driven by congestion processes but 
by large fluctuations of the "normal" traffic behavior. In 
the case of the Internet, the understanding of the physi- 
cal laws governing the nature of traffic is crucial because 
its implications in design, control and speed of the whole 
network 12]. 

In a couple of recent articles, Menezes and Barabasi 
propose a model to understand the origin of fluctuations 
in traffic processes in a number of real world systems, 
including the Internet, the world wide web, and highway 
networks 0, Q|. All of these systems can be repre- 
sented at an abstract level as networks in which pack- 
ets travel from one node to another, packets being real 
data packets or bits in the Internet, files in the world 
wide web, and vehicles in road networks. In particular, 
the authors considered the relationship between the av- 
erage number of packets (fi) processed by nodes during 
a certain time interval, and the standard deviation ai of 
this quantity. They find that there are two classes of 
universality in this relationship for real systems. In the 
Internet, a scales as (/)^^^, whereas a scales as (/) for 
the world wide web and highway networks. Based on 
a stylized model of random walkers throughout the net- 
work, they conclude that this difference is due to the fact 
that the dynamics of the Internet is dominated by "inter- 
nal noise" whereas the dynamics of the world wide web 



users, that is "external noise" . In the abstraction process 
proposed by the authors, they overlook what is probably 
one of the most important factors in the dynamics of traf- 
fic on networks — the limited capacity of nodes to handle 
packets simultaneously, which results in packet-packet in- 
teractions and, eventually, in large fluctuations or even 
network congestion 0, Q . 

In this Letter we show that simple considerations re- 
garding the persistence of packets flowing the network, 
the limitation of nodes to handle information, and the 
time window where statistics are performed, account for 
different scalings of the fluctuations in traflic on complex 
networks. The main results obtained are: (i) Maintain- 
ing the total traffic on the network constant, different 
scaling laws arise depending on the relation between the 
ratio of packet input and the steps these packets per- 
form before disappear; (ii) The time window affects the 
scaling exponent of the fluctuations in such a way that, 
for a small enough time window the scaling trivially sat- 
isfies a ~ ifV^"^ always, no matter the dynamic pro- 
cess. When the time window is large enough, the rest 
of parameters will provide the precise scaling between 
a — 1/2 and a — I where a refers to the scaling exponent 
a ^ if)"', (iii) The effect of the packet-packet interaction 
(queue system) account for different scaling exponents 
as well; (iv) We find that within this framework there is 
not enough evidence for deriving universality classes. We 
have checked the scaling for data flowing on the Abilene 
backbone network, and show that the scaling exponent 
is different from 1/2. 

To understand the origin of the scaling relations for the 
fluctuations in networks, let us consider the behavior of 
a single node — for example, a toll plaza in a highway — 
trying to satisfy demands from users — vehicles arriving 
to the toll-. As we learn from queueing theory [l^j, two 
stochastic processes fully determine the behavior of the 
node: (i) the arrival process by which new packets arrive 
to the node, and (ii) the service process by which the 
node satisfies the demands of the users, that is, forwards 
the packets. The most common queue model corresponds 
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of the packets generation assumes a random (Poisson) 
arrival pattern jig and the service distribution assumes 
a random (exponential) time [l^ . 

Taking into account these considerations, we propose 
to model the traffic process in a complex network of TV 
nodes as N queue systems of type M/M/1, and a random 
walk simulation of packets movement on the network. 
The arrival process of packets to the network is controlled 
by a Poisson distribution with parameter A, each packet 
enters the network at a random selected node. Once the 
packet arrives to the node enters a queue. The delivery 
of the packets in the queue is controlled by an exponen- 
tial distribution of service times with parameter ^. In 
our model, the packets will perform S random steps in 
the network before disappearing, being then S a mea- 
sure of the persistence of packets in the network. This 
dynamics is performed in continuous time, assuming that 
the time expended by packets traveling through a link is 
negligible. 

The system achieves a stationary state whenever the 
arrival rate of packets at each node is lower than or equal 
to the delivery rate, otherwise the system congests. The 
arrival rate at each node i is topology dependent and fol- 
lows a distribution whose mean is X^'^ = BiX where Bi 
is the algorithmic betweenness of node i. Bi is defined 
as the relative number of paths in the network that go 
through node i given a specific routing algorithm . As 
a direct consequence, the node with maximum algorith- 
mic betweenness B^ determines the onset of congestion. 
We will focus on the average number of packets (fi) pro- 
cessed by nodes during a certain time interval P, and the 
standard deviation ai of this quantity. 

Selecting a value of P < = 1/(S*A), we will 

always observe the scaling a ^ {fy^'^, regardless of other 
parameters. Due to the value of P selected, the nodes will 
deliver either one packet or none, at each time interval. 
Suppose that during a number ni of intervals of length 
P the node deliver a packet whereas it docs not deliver 
during a number of intervals uq = rt — rti, where n is the 
number of samples for the statistics. In this situation 
we also have uq 3> ni. Therefore, the average and the 
standard deviation read 

(/)-ni/n (1) 

n ■' 
which can be simplified to 

=[(!-(/))(/)]'/' (2) 

But, in the current scenario, the average flow is (/) <C 
1 and then we recover the a ~ (/)^^^ scaling law. Oth- 
erwise, this argument cannot be applied, and the scaling 
value will be influenced by the rest of parameters of the 
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FIG. 1: (Color online) "Value of the exponent a versus the time 
window of size P, for a fixed Xt-^ = 1/3 and different values 
of the persistence of packets in the network S. The shadowed 
area highlights the region of P in which the exponent a = 
1/2 always appears. The simulation is performed in a scale- 
free network with exponent for the degree distribution 7 = 3 
of 1000 nodes |20l| . "We have observed the same results for 
larger SF networks at a subset of values of P, however the 
computational cost for the whole set of P values used in the 
plot becomes prohibitive. 

In Fig^we show the behavior of the scaling exponent a 
as a function of the time intervals length P over which the 
averages were taken, for a fixed A*''^ = 1/3. We observe 
(shadowed area) that the exponent is always 1/2 when 
the interval length is small enough. Indeed, from the data 
used the exponent 1/2 stands for values of P < 0.01/A*-'^. 

Let us now assume that the sampling of the data is 
performed at intervals of length P 3> l/Xt'^ . In this 
case, we expect the scaling of fluctuations in the system, 
beyond the effect of the sampling process, to be revealed. 
We analyze the behavior of the system varying the rate 
of injection of packets into the system A and the number 
of steps S each node performs before it disappears. We 
first consider that the service rate — > oo. In this case, 
the effect of queues is minimized and then no interaction 
between packets is accounted for. The total traffic T, 
number of packets flowing through the network per unit 
time, is determined by the Poisson process with mean 
(T) = AS*. 

Keeping the total traffic mean (T) fixed, we can con- 
trol the variability of the local traffic incoming to a node 
by varying the values of A and S proportionally. In Fig|21 
we show the scaling exponent transition between a = 1/2 
and a = 1. This plot recovers the results depicted in [l3j, 
although the explanation should be reconsidered in the 
new scenario. The transition of exponent from a = 1/2 
to a = 1 is obtained here simply by increasing the num- 
ber of steps S the packet performs on the network while 
maintaining the mean value of the total traffic (i.e. de- 
creasing proportionally the i njec tion ratio A). This con- 
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FIG. 2: (Color online) Left: Plot a versus (/) for different 
realization of A and S maintaining its product constant. The 
simulation is performed in a scale-free network with exponent 
for the degree distribution 7 = 3 of 1000 nodes |20ll . We fixed 
T = XS = 100. Right: Plot of the a exponent for AS = 100. 
Other values of AS have produced equivalent results, shifted 
to a different region of (/). 



number of steps in the network increases " internal traffic" 
while decreasing the injection of packets is a decrement of 
"external" traffic in this scenario. Nevertheless both re- 
sults are coherent at this point concerning the scaling of 
fluctuations. Our interpretation of this transition is the 
following: for the same total traffic on the network, the 
nature of fluctuations is related to the number of steps 
S the packets perform on the network. When the num- 
ber of steps is small enough the behavior of fluctuations 
is akin a random deposition process independent of the 
topology of the network, X''/ w A. When the number of 
steps in the network grows, the topology induces dynam- 
ical correlations that affect the scaling of fluctuations via 
the algorithmic betweenness, X^^ « XBi. 

We extend the simple model where queues are ne- 
glected, to the more realistic situation when queues are 
persistent. The introduction of queues in the system, in 
our model, is controlled by the parameter /z (rate of ser- 
vice). The possible values of /i are constrained by the 
onset of congestion i.e. /i > A*''^, otherwise congestion 
appears at those nodes with , because of the arrival of 
more packets than those that can be delivered. We inves- 
tigate those values of /z near the onset of congestion to 
reveal the effect of queues in the scaling properties of the 
system. When congestion occurs, the queues correspond- 
ing to those nodes with B^, will have always more packets 
that those than can be delivered in a period P. That 
means that the number of packets delivered by these 
nodes will be controlled exclusively by the service rate /z, 
i.e. the variance scaling with respect to the mean flow at 
these nodes will be again fltted by a = 1/2 corresponding 
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FIG. 3: (Color online) Scaling exponent a as a function of 
the time service for three different time window sizes, and 
for X"/ = 1/3. Shadowed area highlights the region where 
congestion starts at nodes with A*'' = 1/3 . 

of congestion we approach the situation where the scal- 
ing exponent a = 1/2 should be recovered, however the 
possibility that in some periods of time the queues will be 
unoccupied increases as we go away from the congested 
regime, thus a new transition in the scaling exponent as 
a function /i is expected. In Fig. Owe plot the scaling 
exponent transition as a function of /i for a fixed value 
of A^^ = 1/3. In this situation the onset of congestion is 
determined by the critical value fj,c = 1/3. Note that for 
values below /ij. some nodes of the network collapse and 
then gradually the rest of the nodes in the network. In 
this region, shadowed area of Fig. |3| the system enters 
the congestion regime progressively. The transition on 
the scaling exponent depicted in Fig. |3| is also affected 
by the time window size P, we plotted the transition for 
T = 10^, 10^ and lO'*. We observe that as P increases, 
the transition becomes sharper. Indeed in the limit of 
T — > 00 we conjecture that the transition could be dis- 
continuous, and could reflect a first order phase transition 
[21I as observed in other traffic models , although we 
can not claim that this discontinuity will occur sharply 
from 1 to 1/2. 

Up to now, we have show that a simple traffic model 
where the injection of packets to the system follows a 
Poisson distribution, can account for different scaling ex- 
ponents a depending on the parameters A, /i, S and the 
time period P were the statistics are performed. These 
results lead us to suspect that the scaling of fiuctuations 
in real systems must be affected by these parameters as 
well. This cast doubts on the universality predicted in 
[l3| . Indeed, this non- universality has been also claimed 
in the exponent of fluctuations when studying the data 
flow between stocks in NYSE market jl^l- To corrob- 
orate our doubts about universality on the scaling of 
fluctuations in complex networks, we have studied the 
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FIG. 4: Scaling relations between a and (/) for the 112 Abi- 
lene backbone router Interfaces. Analysis performed during 
(a) two days, (b) one week, (c) one month and (d) two months, 
finishing all them in November 15th of 2005. The time win- 
dow size P is fixed to 5 minutes. 



network [2J| that are part of the data also used in 
We collected data from the 112 available router inter- 
faces (links). We gather information of the number of 
packets that exit through each router interface between 
September 15th and November 15th of 2005, at intervals 
of 5 minutes. The scaling a ^ (/)" shows exponents that 
range from a = 0.71 to a = 0.86, significantly different 
from the exponent 1/2 presented in [I3. The interpre- 
tation of these exponents in the context of our stylized 
model is that the Abilene backbone is far from the onset 
of congestion for the interface with maximum algorith- 
mic betweenness, and seems compatible with the mean 
rate of utilization of the interfaces in this backbone that 
is usually below 30%. 

Summarizing, we have presented a simple model of 
traffic in complex networks that capture the essential pa- 
rameters governing the dynamical process. The model 
shows a scaling relationship between a and (/} whose 
exponent depends on the parameters considered as well 
as on the time window in which the statistics are per- 
formed. Moreover we have shown that the corresponding 
exponent for the scaling of fluctuations in the Internet 
Abilene backbone network is different from 1/2 as stated 
in previous works, corroborating by exclusion that the 
universality on the scaling of fluctuations in complex net- 
works should be questioned. 
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